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INTRODUCTION 


This  report  describes  the  development  of  a  self-tuning  control 
algorithm  suitable  for  application  to  many  processes,  including  those  in 
heating,  ventilating,  and  air  conditioning  (HVAC)  systems.  The  controller 
was  developed  to  reduce  the  effort  required  to  set  up  and  maintain 
controls  for  HVAC  systems  and  to  assure  optimum  control  (and,  consequently, 
maximum  operating  efficiency)  over  a  wide  range  of  conditions. 


BACKGROUND 

Small  independently  operating  digital  control  systems  were  intro¬ 
duced  by  the  building  controls  industry  a  few  years  ago  as  an  alternative 
to  the  traditional  pneumatic  and  analog  electronic  control  systems  for 
control  of  heating,  ventilating,  and  air  conditioning  equipment.  These 
new  control  systems,  called  direct  digital  control  (DDC)  systems,  are 
rapidly  gaining  acceptance  by  building  designers,  owners,  and  operators 
because  of  their  accuracy,  versatility,  reliability,  and  cost  effective¬ 
ness.  The  benefits  of  DDC  may,  however,  come  at  the  expense  of  increased 
complexity  in  adjusting  the  control  system.  For  example,  most  pneumatic 
controllers  provide  only  proportional  control,  while  a  DDC  system  can 
readily  provide  the  more  accurate  proportional-integral-derivative  (PID) 
control.  However,  to  gain  the  benefits  of  PID  control  requires  consid¬ 
erably  more  effort  to  determine  the  correct  values  of  control  parameters, 
such  as  gains,  to  input  into  the  controller.  In  addition,  the  optimum 
control  parameters  for  one  mode  of  operation,  such  as  heating,  may  not 
be  the  optimum  parameters  for  another  mode  of  operation,  such  as  cooling. 
If  the  control  parameters  for  a  DDC  system  are  not  chosen  with  care,  the 
resulting  control  may  well  be  worse  than  that  obtainable  with  the  simplest 
pneumatic  control  system.  Therefore,  it  is  very  desirable  to  develop  a 


controller  methodology,  or  algorithm,  which  could  automatically  determine 
what  the  optimum  control  parameters  should  be  and  change  the  control 
parameter  values  as  operating  conditions  require.  A  few  years  ago,  the 
computational  requirements  of  such  a  self-tuning,  or  adaptive,  controller 
would  have  required  a  sizable  computer  to  implement.  Today,  however, 
the  computational  requirements  of  the  algorithm  presented  below  are 
within  the  capabilities  of  many  small  DDC  systems. 

DISCUSSION 

Approach 

The  guidelines  for  developing  the  controller  methodology  were: 

(1)  the  method  must  not  require  the  input  of  a  detailed  mathematical 
model  of  the  process  being  controlled,  and  (2)  the  method  should  not  be 
sensitive  to  variations  in  time  delays.  Although  adaptive  control 
algorithms  have  been  developed  and  applied  to  the  control  of  HVAC  compo¬ 
nents  (e.g.,  Park  and  David,  1982),  these  previous  controller  designs 
required  input  of  a  model  of  the  process  being  controlled.  The  authors 
felt  that  this  modeling  requirement  would  severely  restrict  the  appli¬ 
cations  of  a  self-tuning  controller  because  it  would  be  difficult  for 
control  engineers  to  develop  the  necessary  models  for  each  controller 
installation.  Also,  it  was  desired  that  the  self-tuning  algorithm  be 
insensitive  to  variations  in  system  time  delays,  such  as  response  lag, 
as  these  delays  would  also  be  different  at  each  controller  installation 
and  change  with  the  process  dynamics. 

The  first  step  in  the  development  of  a  self-tuning  controller  for 
HVAC  applications  was  to  conduct  a  thorough  review  of  past  work  on 
adaptive  control  and  select  the  most  promising  approach  to  further 
analysis.  This  work  was  reported  by  Seborg  et  al.  (1983). 

Based  on  the  guidelines  and  review  of  past  work,  a  self-tuning 
controller  was  designed  by  I-Lung  Chien  (Chien,  1985).  This  controller 
is  described  below  in  detail.  The  controller  was  then  evaluated  using  a 
computer  simulation  of  a  hot  water  coil.  The  performance  of  the  self-tuning 


controller  was  compared  to  that  of  an  optimally  tuned  proportional-integral 
(PI)  controller  operating  under  the  same  conditions.  The  optimum  controller 
settings  for  the  PI  controller  were  determined  using  analytical  and 
experimental  methods  and  were  manually  input  into  the  controller.  On 
the  basis  of  the  encouraging  results  of  the  computer  simulations,  experi¬ 
ments  on  an  actual  hot  water  coil  system  were  conducted. 

Controller  Design 

The  self-tuning  controller  is  based  on  the  following  design  strategy. 
Parameters  in  a  simple  dynamic  model  of  the  process  are  updated  on-line. 

The  controller  setting  calculations  are  then  based  on  the  current  parameter 
estimates.  The  dynamic  model  is  a  simple,  empirical  model  rather  than  a 
detailed  physical  model.  Thus  the  user  need  only  specify  the  model 
order  and  the  estimated  time  delay.  Typically,  a  second-order  plus  time 
delay  model  is  adequate.  The  model  parameters  are  updated  at  each 
sampling  instant  from  input-output  data  using  standard  least  squares 
techniques.  The  controller  settings  are  then  calculated  in  a  straight¬ 
forward  fashion  from  the  model  parameters  so  as  to  minimize  a  quadratic 
cost  function. 

The  starting  point  of  the  analysis  is  a  single-input/single-output 
(SISO) ,  discrete-time  process  model  (see  List  of  Symbols  and  Figure  1 
for  notation). 
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where  the  polynomials  A  and  B  are  expressed  in  terms  of  the  backward 
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and  the  time  delay  represented  by  k  is  expressed  as  a  multiple  of  the 
sampling  period  (k  >  1),  while  t  denotes  the  sampling  instant,  t  =  0,  1, 
2,  ....  It  is  also  assumed  that  the  poles  of  A(z  *)  lie  inside  the  unit 
circle  in  the  complex  z-plane,  i.e.,  that  the  system  is  open-loop  stable 
Consider  the  Smith  predictor  approach  to  time  delay  compensation 
(Smith,  1957).  Figure  1  shows  the  block  diagram  of  the  discrete-time 
Smith  predictor.  Transfer  function  G' (z  ^)  denotes  the  process  model 
without  the  time  delay,  i.e.,  with  k  =  1: 
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Adapting  the  Smith  predictor  to  accommodate  parameter  variations 

requires  estimating  the  a^  and  b^  parameters  of  the  process  model  as 

well  as  the  time  delay,  k.  However,  as  discussed  by  Kurz  and  Goedecke 

(1981)  and  Vogel  (1982),  time  delays  are  difficult  to  estimate  on-line 

directly,  and  thus  this  approach  should  be  avoided.  In  an  alternative 

approach,  Vogel  and  Edgar  (1982)  assumed  that  the  minimum  and  maximum 

values  of  the  expected  time  delay,  k  .  and  k  ,  are  available.  They 

min  max 

then  performed  the  parameter  estimation  by  rewriting  Equation  1  with  a 
sufficient  number  of  terms  in  the  numerator  of  the  process  transfer 
function  model  to  include  implicitly  any  extra  time  delay  in  excess  of 
k  .  Using  their  approach,  the  process  model  becomes 
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Note  that  the  variation  in  the  process  time  delay  is  accounted  for 

implicitly  when  leading  coefficients  of  B  are  zero,  since  B  can  be 
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expressed  as 
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Next,  the  process  model  without  time  delay,  G'(z  ),  is  redefined 
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G^(z  1)  is  chosen  in  this  fashion  because  it  eliminates  the  time  delay 
in  the  process  model  without  requiring  explicit  knowledge  of  k,  and  yet 
it  contains  the  same  denominator  as  G  (z  *)  and  has  the  same  steady-state 

P  -1 

gain.  It  is  important  to  note  that,  with  the  above  choices  of  G  (z  ) 

-1  P 

and  G^(z  ),  the  difficult  problem  of  explicit  on-line  time  delay  estima¬ 
tion  is  completely  avoided. 

Note  also  that  G^(z  1)  provides  only  an  approximate  prediction  of 
the  undelayed  process  output.  Vogel  (1982)  shows  that  this  time  delay 
compensator,  even  with  the  inherent  approximation,  generally  has  a  more 
beneficial  effect  on  the  closed-loop  performance  and  stability  than  the 
standard  Smith  predictor. 

In  principle,  the  design  of  the  controller  G  (z  *)  in  Figure  1  can 

-1  c 

be  based  on  the  undelayed  process  model  G^(z  )  since  the  time  delay  is 

eliminated  from  the  characteristic  equation.  The  SISO  self-tuning 

controller  of  Clarke  and  Gawthrop  (1975,  1979)  for  G' (z  *)  can  be 

P 

obtained  as: 
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where  y(t)  is  the  measured  output,  u(t)  is  the  manipulated  input,  w(t) 
is  the  setpoint,  and  d  is  a  constant,  steady-state  bias  response  for  a 
zero  input  signal.  E  and  F  are  polynomials  in  z  *  which  satisfy  the 
identity, 


P  =  EA  +  z“k  F 

and 
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For  the  undelayed  process  model,  k  =  1  and  thus  E  =  1.  P,  Q,  and  R  are 
user-specified  polynomials  to  achieve  a  desired  response  time  and  degree 
of  oscillation  for  the  closed-loop  system.  Guidelines  for  the  choice  of 
these  design  parameters  can  be  found  in  Chien  (1985). 

Figure  2  shows  the  block  diagram  when  the  self-tuning  controller  is 
used  in  conjunction  with  Smith-predictor-type  time  delay  compensation 
(STC-TDC).  The  resulting  control  law  can  be  expressed  as 


u 
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p  \  +  Q  A  -  F  be 


min 


where 


e  =  Rw  -  Fy  -  d 


(15) 


(16) 


In  summary,  at  each  sampling  instant,  the  following  calculations 
are  performed  for  the  new  STC-TDC: 


1.  Read  new  values  for  the  output  y(t)  and  setpoint  w(t). 


2.  Update  model  parameters  A,  B  ,  and  d  by  using  a  parameter 

L 

estimation  scheme  such  as  the  recursive  least-squares  algorithm. 

3.  Obtain  F  from  Equation  12. 

4.  Calculate  the  new  control  action  from  Equation  15. 

Chien  (1985)  also  describes  a  load  disturbance  rejection  scheme  to 
improve  the  load  response  during  upset  periods. 

Simulation  Results 


A  heating,  ventilating,  and  air  conditioning  process  model  developed 
by  Todtli  et  al.  (1982)  has  been  used  to  evaluate  the  performance  of 
this  STC-TDC.  This  process  consists  of  an  air  heater  operating  under 
conditions  of  varying  air  inlet  temperature,  air  flow  rate,  and  water 
inlet  temperature.  A  schematic  diagram  of  the  air  heater  is  shown  in 
Figure  3.  The  variables  in  Figure  3  are  defined  as: 


u  =  normalized  valve  position  (0  <  u  <  1) 

T^  =  water  inlet  temperature,  °C 

T' .  =  water  inlet  temperature  at  the  entrance  of  the  air 
W1  heater,  °C 

T  =  water  outlet  temperature,  °C 

WO 

T  =  air  inlet  temperature,  °C 

Tao  =  air  outlet  temperature,  °C 

m  =  mass  flow  rate  of  water  after  the  control  valve,  kg/sec 
w 

m’  =  mass  flow  rate  of  water  at  the  entrance  of  the  air  heater, 
W  kg/sec 

m^  =  mass  flow  rate  of  air,  kg/sec 
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From  a  control  point  of  view,  the  process  has  four  input  variables:  u, 

T  . ,  m  ,  and  T  . ;  u  is  the  single  manipulated  variable  and  the  other 
ai  a  wi 

three  variables  are  considered  to  be  disturbances.  The  output  variable 
that  has  to  be  controlled  is  T 

ao 

Based  on  experimental  data  for  a  particular  commercial  air  heater, 

the  process  model  that  describes  its  static  and  dynamic  behavior  was 

developed  by  Todtli  et  al.  (1982).  A  detailed  derivation  of  the  process 

model  can  be  found  in  Chien  (1985). 

The  nominal  operating  point  is  selected  to  be:  m  =  0.388  kg/sec, 

a 

T^^  =  90°C,  T  =  0°C,  and  u  =  0.262.  The  resulting  air  outlet  tempera¬ 
ture  at  steady  state  is  20°C.  Figure  4  shows  the  open-loop  response  for 
changes  in  u  at  this  operating  point,  and  Table  1  shows  the  time  con¬ 
stants  and  gains  for  these  changes  in  u. 

Figure  5  shows  the  operating  window  for  this  HVAC  process  in  the 

operating  regime:  -20°C  §  T  .  §  20°C,  0.2  kg/sec  §  m  §1.0  kg/sec,  and 

ai  a 

T  .  =  90°C. 
wi 

The  HVAC  process  model  described  above  was  simulated  to  examine  the 
performance  of  the  STC-TDC  in  comparison  with  a  PI  controller.  For 
control  purposes,  a  sampling  period  of  10  seconds  was  chosen;  process 
noise  was  simulated  by  adding  a  Gaussian  noise  signal  with  a  standard 
deviation  of  0.5°C  to  the  air  inlet  temperature,  T  .  Setpoint  changes 
were  made  at  t  =  300  seconds  and  t  =  1,300  seconds,  while  a  load  change 
was  made  at  t  =  800  seconds. 

A  linearized  model  of  the  above  HVAC  process  valid  around  the 
nominal  operating  point  was  used  to  determine  Ziegler-Nichols  settings 
for  the  PI  controller.  The  resulting  PI  controller  settings  were  cal¬ 
culated  to  be  K  =  0.0793  1/°C  and  tT  =  83.2  seconds.  For  STC-TDC,  the 
P  I 

assumed  process  model  used  had  n.  =  1 ,  r  =  3 ,  n„  =  0,  and  k  .  =1. 

A  C  ^pin 

Also,  the  design  parameters  were  chosen  to  be  P  =  1  -  0.7z  ,  q  =  1  -  z  , 

and  R  =  0.3. 

Figures  6  and  7  show  the  performance  of  the  PI  and  STC-TDC  controller 
when  the  air  inlet  temperature  increases  by  5°C.  The  PI  controller 
becomes  oscillatory  after  the  load  disturbance,  while  the  STC-TDC  per¬ 
forms  very  well  for  both  setpoint  and  load  changes.  Figures  8  and  9 
show  the  performance  of  the  two  controllers  when  air  mass  flow  rate 
provides  a  load  disturbance.  It  is  assumed  that  the  process  time  delay 
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decreases  to  10  seconds  when  air  mass  flow  rate  increases  20%.  When  the 
time  delay  decreases,  the  PI  controller  yields  about  the  same  response 
as  the  STC-TDC,  as  shown  in  Figures  8  and  9.  However,  Figures  10  and  11 
show  that  when  time  delay  increases,  the  PI  controller  becomes  unstable 
after  the  load  change  occurs,  while  the  STC-TDC  still  gives  perfect 
response . 


EXPERIMENTAL  RESULTS 

Apparatus 

The  experimental  apparatus  is  diagrammed  in  Figure  12.  A  portable 
computer  was  purchased  to  serve  as  both  a  programmable  controller  and  as 
a  data  logging  device  (Figure  13).  The  computer  was  equipped  with  a 
high  performance  analog  and  digital  input/output  (AD/DA)  board  to  enable 
data  to  be  read  into  the  computer  and  enable  the  output  of  the  control 
program  to  affect  the  process.  Experimental  results  were  recorded  on 
one  of  the  computer's  disk  drive  units  to  facilitate  later  analysis. 
Results  were  also  printed  to  provide  a  hardcopy  record  of  the  experiments 
A  real  time  graphics  program  was  written  to  display  the  results  of  the 
experiment  while  it  was  in  progress. 

The  temperature  sensors  were  1,000-ohm  platinum  resistance  tempera¬ 
ture  devices  (RTD)  mounted  on  the  ends  of  duct  probes.  Current  loop 
transmitters  located  near  the  sensors  converted  the  temperature-dependent 
voltage  signals  from  the  RTDs  to  electrical  currents  to  reduce  the 
problems  caused  by  electrical  noise  in  the  input  signals.  The  current 
loop  transmitters  were  connected  to  precision  load  resistors  on  the 
AD/DA  board  by  shielded  cable. 

The  output  voltage  from  the  AD/DA  board  was  converted  into  an  air 
pressure  signal  by  means  of  an  electrical-to-pneumatic  transducer.  The 
air  pressure  signal  operated  the  hot  water  control  valve. 


Method 


A  computer  program  was  written  to  experimentally  verify  the  perfor¬ 
mance  of  the  self- tuning,  time-delay-compensating  (STC-TDC)  controller 
algorithm.  The  computer  program  actually  was  capable  of  performing 
three  types  of  controller  experiments:  open-loop  tests,  closed-loop  PID 
tests,  and  the  self-tuning  controller  tests.  Open-loop  tests  were 
conducted  on  the  experimental  apparatus  to  determine  the  response  charac¬ 
teristics  of  the  hot  water  valve/coil  system.  These  characteristics 

were  used  to  determine  the  optimum  values  of  proportional  gain  and  reset 

rate  for  the  PI  controller  tests.  The  hot  water  coil  was  controlled  by 

an  optimized  PI  controller  and  the  results  were  used  as  a  basis  for 

judging  the  performance  of  the  self-tuning  controller.  Experiments  were 
then  conducted  using  the  self-tuning  controller.  Experiments  were 
performed  for  two  cases:  (1)  the  load  on  the  coil  (as  measured  by  the 
coil  inlet  temperature)  is  constant  but  the  setpoint  (as  measured  by  the 
coil  discharge  temperature)  is  changed,  and  (2)  both  the  load  on  the 
coil  and  the  setpoint  are  changed. 

Results 


Eleven  separate  experiments  were  conducted  using  an  air  handler  at 
the  HVAC  Test  Facility  (Figure  14)  at  the  U.S.  Army  Corps  of  Engineers 
Construction  Engineering  Research  Laboratory  (CERL)  in  Champaign,  Illinois 
Occasional  noise  in  the  input  data  for  the  first  few  experiments  had  a 
serious  detrimental  effect  on  the  performance  of  the  STC-TDC  algorithm. 
After  the  usual  fixes  for  problems  of  signal  noise  (i.e.,  connecting 
common  circuit  grounds,  adding  by-pass  capacitors)  failed  to  eliminate 
the  noise  spike  problem,  the  computer  program  was  modified  to  add  a 
digital  filter  to  the  input  data.  The  digital  filter  cured  the  noise 
problem.  The  second  problem  with  the  experimental  test  activity  was 
that  the  valve  used  to  control  the  flow  of  hot  water  to  the  heating  coil 
was  greatly  oversized.  As  a  consequence,  the  valve  operated  at  only 
about  10  to  15%  of  capacity  at  the  design  point,  rather  than  the  standard 
50%  of  capacity  at  the  design  point.  The  resultant  reduced  throttling 
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range  made  it  difficult  to  control  the  coil  with  any  precision.  For 
reasonable  values  of  changes  in  setpoint,  the  valve  was  either  open  or 
closed,  with  the  result  that  any  control  method  behaved  like  on-off 
control.  The  problem  of  the  oversized  valve  was  minimized  by  conducting 
experiments  using  only  comparatively  small  changes  in  setpoint.  The 
comparatively  small  temperature  range  over  which  the  coil  discharge 
temperature  was  controlled  emphasized  the  effect  of  another  variable: 
hot  water  supply  temperature.  Hot  water  was  supplied  to  the  coil  by  a 
steam  convertor  and  the  temperature  of  the  water  was  controlled  by 
on-off  regulation  of  the  flow  of  steam  to  the  convertor.  The  temperature 
of  the  supply  water  fluctuated  by  about  0.5°C  which,  unfortunately,  was 
about  the  magnitude  of  the  changes  in  setpoint. 

Figures  15  and  16  present  the  results  of  an  experiment  to  determine 
the  performance  of  a  PI  controller  under  the  conditions  of  setpoint 
change  but  no  load  disturbance.  Figure  15  shows  controller  performance; 
Figure  16  shows  coil  inlet  temperature.  The  apparent  oscillation  in 
temperature  about  the  setpoint  is  attributed  to  the  fluctuating  tempera¬ 
ture  of  the  supply  water. 

Figures  17  and  18  show  the  results  for  a  PI  controller  experiment 
in  which  the  coil  inlet  air  temperature  was  suddenly  changed. 

Figures  19  through  22  present  results  for  the  STC-TDC  controller 
comparable  to  those  presented  in  Figures  15  through  18  for  a  PI  con¬ 
troller.  Figure  19  illustrates  how  the  STC-TDC  controller  "learns"  as 
information  about  system  response  to  disturbances  is  received  and  pro¬ 
cessed.  Notice  (in  Figure  19)  how  the  response  of  the  STC-TDC  controller 
to  the  first  setpoint  change  (at  time  equals  10  minutes)  is  quite  poor, 
the  controller  response  to  the  second  setpoint  change  (at  time  equals 
30  minutes)  is  better,  and  the  controller  response  to  subsequent  setpoint 
changes  (at  time  equals  60  minutes)  is  comparable  to  that  of  an  optimally 
tuned  PI  controller  (see  Figure  15). 

Figure  21  presents  the  results  of  an  experiment  to  evaluate  the 
response  of  the  STC-TDC  controller  to  a  disturbance  in  the  coil  load. 

At  time  equals  60  minutes,  the  outside  air  dampers  were  opened  to  admit 
cold  air  and  lower  the  coil  inlet  temperature.  As  Figure  21  shows,  the 
response  of  the  STC-TDC  was  less  than  satisfactory.  The  STC-TDC  provided 
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good  initial  recovery  from  load  upset  (using  the  load  rejection  scheme 
described  by  Chien  (1985)).  However,  the  STC-TDC  tended  to  switch  back 
to  the  standard  STC-TDC  control  law  too  soon,  which  introduced  an  addi¬ 
tional  upset.  It  is  believed  that  this  problem  can  be  eliminated  by  a 
simple  change  to  the  computer  code.  However,  additional  experimental 
tests  are  necessary  to  confirm  this  hypothesis  (see  Chapter  4  of  Chien 
(1985)). 

Before  the  experiment  could  be  repeated  with  a  modified  STC-TDC, 
the  hot  water  convertor  at  the  HVAC  Test  Facility  broke  down  and  the 
remainder  of  the  experiment  had  to  be  postponed  until  the  hot  water 
convertor  was  repaired. 

Funding  for  the  work  described  in  this  report  expired  before  addi¬ 
tional  experimental  work  could  be  performed. 

It  should  be  noted  that  the  STC-TDC  has  been  successfully  applied 
to  a  pilot-scale  multicomponent  distillation  column  by  Chien  (1985). 

CONCLUSIONS 

An  algorithm  for  a  self-tuning  controller  was  developed  and  tested. 
The  control  methodology  has  potential  for  wide-spread  application  in 
industrial  processes  such  as  control  of  the  apparatus  of  heating,  venti¬ 
lating,  and  air  conditioning  systems. 

The  controller  was  developed  to  reduce  the  effort  required  to  set 
up  and  maintain  the  modern,  more  complex,  digital  control  systems  being 
installed  on  HVAC  equipment.  A  controller  that  could  adjust  itself  to 
maintain  optimum  HVAC  operating  conditions  would  result  in  improved 
plant  efficiency  and,  possibly,  improved  equipment  reliability.  Mainte¬ 
nance  expenditures  for  controller  adjustments  would  be  eliminated. 

Although  the  controller  has  been  demonstrated  to  work  well  in  some 
applications,  additional  experimental  work  is  required  to  verify  its 
applicability  to  the  variety  of  controlled  devices  found  in  HVAC  systems 


RECOMMENDATIONS 


It  is  recommended  that  additional  experimental  work  be  performed  to 
evaluate  the  performance  of  the  STC-TDC  algorithm  for  control  of  a  hot 
water  coil  subject  to  changes  in  coil  inlet  air  temperature.  The  experi 
mental  evaluations  performed  using  a  hot  water  coil  should  be  repeated 
using  a  chilled  water  coil.  It  is  also  recommended  that  the  STC-TDC 
controller  algorithm  be  evaluated  for  control  of  a  two-phase  coil,  such 
as  a  steam  coil  or  a  direct  expansion  coil,  and  a  variable  speed  fan. 

If  these  tests  prove  successful,  the  algorithm  should  be  installed  on  a 
commercially  available  model  of  a  digital  controller  for  extended  field 
evaluation. 
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Figure  1.  Block  diagram  of  discrete-time  Smith  predictor 


Figure  2.  Block  diagram  of  the  STC  with  time  delay  compensation 
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10.  PI  controller  for  setpoint  changes  when  air  mass  flow  rate  (m  ) 
decreases  20%  and  time  delay  increases  to  30  sec  at  t  =  800  sec 


Diagram 


Figure  13-  Portable  computer  used  as  programmable  controller  and  data  logger. 
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Figure  15.  Coil  discharge  temperature  data  from  experiment  PI.l 


Figure  16.  Coil  inlet  temperature  data  from  experiment  PI.l 
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